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Abstract 

The interaction of electric field with charged domain walls in ferroelectrics is theoretically ad- 
dressed. A general expression for the force acting per unit area of a charged domain wall carrying 
free charge is derived. It is shown that, in proper ferroelectrics, the free charge carried by the wall 
is dependent on the size of the adjacent domains. As a result, it was found that the mobility of such 
domain wall (with respect to the applied field) is sensitive to the parameters of the domain pattern 
containing this wall. The problem of the force acting on a planar charged 180° domain wall normal 
to the polarization direction in a periodic domain pattern in a proper ferroelectric is analytically 
solved in terms of Landau theory. It is shown that, in small applied fields (in the linear regime), 
the forces acting on walls in such pattern increase with decreasing the wall spacing, the direction of 
the forces coinciding with those for the case of the corresponding neutral walls. At the same time, 
for large enough wall spacings and large enough fields, these forces can be of the opposite sign. It 
is shown that the domain pattern considered is unstable in a defect-free ferroelectric. The poling 
of a crystal containing such pattern, stabilized by the pinning pressure, is also considered. It is 
shown that, except for a special situation, the presence of charge domain walls can make poling 
more difficult. It is demonstrated that the results obtained are also applicable to zig-zag walls 
under the condition that the zig-zag amplitude is much smaller than the domain wall spacing. 
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I. INTRODUCTION 



Typically, at domains walls in ferroelectrics, the normal component of the electrical dis- 
placement is conserved with high accuracy This is the so-called condition of electrical 
compatibility 1 . When a ferroelectric behaves as a perfect insulator, any appreciable viola- 
tion of this condition will imply the appearance of bound charge at the wall and a macro- 
scopic electric field in the adjacent domains. Such a field, at best, can strongly increase 
the energy of the system, but very often this field is expected to be strong enough to fully 
suppress the ferroelectricity in the sample. At the same time, the walls with an essential 
violation of the electrical compatibility - the so-called charged domain walls - have been 
observed in different materials such as lead titanate PbTiOar - -, PbfZr^Tii^jOs (PZT)^, 
lead germanate PbsGeaOn (PGO)->£, and lithium niobate^. It is believed that such walls 
can exist due to the compensation of its bound charge by free electron^"- or ion charges^. 
Here, the situation with the electron compensation is of special interest. Specifically, as 
was recognized many yeas ago, the concentration of free carriers in properly compensated 
domain walls can readily achieve the metallic level^^. Thus, one can treat such a wall as 
a highly conductive interface. In principle, one can tune its conductivity by controlling the 
wall orientation^ and change its position in the crystal by the application of a dc electric 
field. In this context, mobility of charged domain walls is an important issue. Additionally, 
this issue deserves attention in view of poling samples containing charged domains walls. 

The dynamics of a charged domain wall has not been completely understood and theo- 
retical results are limited. Here, one can mention only a classical paper by Landauer— and 
a recent paper by Mokry et al.— Landauer indicated that the compensation of the bound 
charge on a wall will lead to a reduction of the pressure acting on it in the presence of a dc 
electric field, whereas Mokry et al. qualitatively describled this effect in terms of Landau 
theory. A closer analysis of the problem, however, shows that it misses a more involved 
theoretical treatment. For example, a straightforward application of the results by Mokry 
et al.— to the case of a "zig-zag" charged wall leads to a paradox. 

The goal of this paper is to revisit the problem of mobility of charged domain walls. The 
paper is organized as follows. 

In Sect. [Til we discuss and solve the aforementioned paradox. Here, we obtain the general 
expression for the local force acting on the domain wall. It is shown that not only the normal 
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component of the force (treated in Ref.— ), but also the tangential one should be considered. 
In Sect. II III we considered the force acting of a planar charged 180° domain wall normal to 
the polarization direction in a periodic domain pattern in a proper ferroelectric. We have 
shown that, for a fixed applied field, the force acting on the domain wall is a function of 
the structure period. The results obtained are applied to the analysis of poling of such a 
structure. Section [TV] is devoted to a discussion of the implications of the results obtained 
to the experimental situation and available experimental data. In Appendix, supporting 
calculations for Sect. |TTT]are presented. 



II. LOCAL FORCE DENSITY ON A FERROELECTRIC DOMAIN WALL 

In this section we will consider the problem of the force acting on a charged domain 
wall placed in an electric field. Originally, this problem was addressed by Mokry et al.— , 
who obtained a relation linking the pressure acting on the field (the surface density of the 
normal component of the force) with values of polarization and electric field in the adjacent 
domains and surface density of the free charge carried by the walls. However, one finds 
a situation where a straightforward application of this result leads to a paradox. This 
can be demonstrated for the case of the force acting on charged "zig-zag" walls, which are 
typically observed in experiments^"-^. Consider a fully compensated (the total (bound+free) 
charge on the wall equal to zero) zig-zag wall in a parallel plate capacitor of a thickness h, 
shown in Fig. [TJ Let us do this in hard ferroelectric approximation where the electrical 
displacement inside the domains is presented as a sum of the spontaneous polarization and 
the linear dielectric response to the electric field with permittivity £f . The bottom electrode 
is grounded and the top electrode is connected to a source of constant voltage —U. In this 
case, the pressure on each flat element of the wall was obtained by Mokry et al.— as: 

p = 2P E cos 2 6, (1) 

where Pq is the spontaneous polarization, Eq = U/h is the applied field, 9 is the angle 
between the wall and the applied electric field. The directions of the corresponding forces 
are shown in Figfjji. As is clear from this figure, such expression for the pressure leads to a 
nonzero total force acting on the zig-zag structure. 

On the other hand, we can obtain this force directly from the principle of virtual displace- 
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FIG. 1. Zig-zag wall in a parallel plate capacitor. The pressure acting on each segment leads to 
non-zero force acting on the wall (a). During the virtual displacement of the "zig-zag" domain 
wall the charge Q on the top electrode remains unchanged. Thus the force on the wall, calculated 
directly with the principle of virtual displacement, is zero (b). 

ment. This method leads to a zero force acting on the fully compensated zig-zag domain 
wall. 

Indeed, according to this principle, the work 5W done by the external force F is equal 
to the variation of the proper thermodynamic function SG during virtual displacement 8u: 
5W = F5u = 5G. For the system considered, in the hard ferroelectric approximation, this 
thermodynamic function can be denned as^: 



where Q is the charge on the top electrode, E is the electric field. In view of full charge 
compensation, the electric field inside the ferroelectric is equal to U / h and does not change 
with virtual displacement. The virtual displacement of the "zig-zag" wall will not affect 
the charges on the electrodes (see Fig. QJb)), thus the thermodynamic function G will not 
change with the virtual displacement. This implies that the force acting on the wall is zero. 

The above paradox can be resolved if we take into account the fact that the local force 
is not necessarily perpendicular to the domain wall. Thus to calculate the resultant force, 
acting on an element of the domain wall, it is not sufficient to know the pressure. The 
component of the force parallel to the domain wall should also be found. As we will see 
later, in the case of the fully compensated wall, the consideration of all components of the 
force density will give the resultant force equal to zero. 

Let us obtain the general expression for the force density on the domain wall. We will 
start from scratch using the principle of virtual displacement. We can not use the formula 




(2) 
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for the generalized stress tensor from textbook by Landau and Lifshitz 15 , as was done in 
Ref.— . In fact, no derivation of this formula is available. In addition, one can show that it 
leads to results which are in contradiction with the results obtained with the principle of 
virtual displacement. 

In our analysis we will use the same approximations as those used by Mokry et al.— . First, 
the internal structure of domain wall is neglected. Second, the free charge is moving together 
with the domain wall. Third, we will not consider the forces arising due to an imbalance of 
the elastic energy, which is a suitable approximation for non-ferroelastic domain walls. 

In order to obtain a general formula for the local force density on domain wall, we follow 
the principle of virtual displacements, as applied to the domain wall of arbitrary shape. 
We will consider the problem at fixed potentials on the conductors ip^} . We will use the 
thermodynamic function G of the general form: 

dv-J2^Q ( l (3) 

i 

where the first term represents the part of the free energy, $(P), associated with the fer- 
roelectric part of polarization P and the electric field energy integrated over the volume V 
of the ferroelectric (E is the electric field); the second term represents the subtracted work 
of the electric sources (Q^) are the charges on the conductors). Here, is the background 
permittivity, which was for simplicity taken as isotropic. Thus, in our notation, the vector 
of electric displacement reads: 

A = e h Ei + Pi. (4) 

More or less straightforward analysis in terms of virtual displacement, given in Appendix |A~1 
leads to the following expression for the force acting on a unit area of the the domain wall: 

f k = [[$(P)]] n k - % [[Pi]] n k + E k a u (5) 

where [[Z]] = Z® - Z« denotes the jump of the quantity Z on the domain wall, E = 
(^K 1 ) _|_ E^)/2, and n k is a unit vector, normal to the domain wall, directed from domain 
(1) to domain (2). The upper suffix is used to designate the domains. 

First two terms in Eq. (pjj) are the same as for the neutral domain wall. The third term 
is the Coulomb force acting on the free charge on the domain wall. 

5 



G 



v 



$(P) + l -e h E> 



We can rewrite Eq. (J5J) as the sum of the components perpendicular and parallel to the 
domain wall as: 

f k = {[MP)}) - % ([[P.]] - a fni )^n k + \E k - E 3 n 3 n k } a { . (6) 

The first term in the curly brackets is the pressure acting on the domain wall, which is the 
same as that found before by Mokry et al.— : 

V = [[®(P)]\ ~ % ([[Pi\] - a f m) . (7) 

In this particular case the results, obtained with the principle of virtual displacement and 
with the generalized stress tensor from the Landau book— are the same^. 
The second term in Eq. represents the force parallel to the domain wall 

/|| = £||*f, (8) 

where En k = E k — Ejrijn k is the component of the average electric field parallel to the 
domain wall. The force /|| is the tangential component of the Coulomb force acting on the 
free charges on the wall. 

Let us use the result obtained to find the force acting on a fully compensated zig-zag 
domain wall, the situation where we have come across a paradox. As before, we consider the 
problem in the hard ferroelectric approximation. We consider the model of a thin capacitor, 
with the thickness h, much smaller than the capacitor width. It is the same model as used 
to obtain the well-known formula for the capacitance of the thin parallel-plate capacitor. 
In this model, to exclude the edge effects and to consider a uniform electric field inside 
the capacitor, we should deal with the force per unit area of the electrode in the limit of 
infinitely wide capacitor instead of the total force acting on the wall. 

The free charge density on the completely compensated " head-to-head" domain wall, that 
is equal to the jump of the normal component of the spontaneous polarization, reads: 

a f = -2P sin 9. (9) 

Thus, from Eq. ([6]) one finds the force acting on a flat segment of the domain wall as 

F k ={2P E cos 2 9 n k - 2P sin 9{E ok - E sin 9 n k )}S = 

( E 0k \ ( 10 ) 

--2P E S [ n k - sm6»— 



where S is the area of the wall segment. This force is perpendicular to the direction of the 
electric field, i.e. it is parallel to the capacitor electrodes. Indeed, one can check that 



F k E ok = 2P E S (n k E ok - E sin 6) = 0. (11) 



The second term in brackets in Eq. ( flOl) is parallel to the electric field, thus it is not 
contribute to the component of the force parallel to the capacitor electrodes, and the force 
can be found as 

F k = 2P E Sd k , (12) 

where d k is the projection of n k on the plane of the electrodes. 

Thus the absolute value of the force can be found from Eq. ffl2|) as: 

|F| = 2P E Scos6. (13) 

This force is acting in the direction of a vector d k , i.e. from domain 1 to domain 2. Thus 
forces acting on the neighbor segments are oppositely directed (see Fig. [2]). Thus the total 
force acting on the zig-zag wall containing n segments can be found as the sum the the forces 
acting on each segment 

n n 

F = 2P E (-^)% cosOi = 2P E l £ (z t - z^) = 2P E l(z n - z ), (14) 

i=l i=0 

where the z-axis is directed perpendicular to the electrodes, Zi is the coordinate of the zig- 
zag vertex, / is the capacitor width in the direction, perpendicular to the figure plane (Fig. 
[2]). The force obtained depends only on the coordinates of the vertices on the edges of the 
capacitor. Considering the density of this force per unit area of the electrode we get 

F 2P E l{z n - zq) ^ 2P E h 

Del IW W 

where w is the width of the capacitor (see Fig. [2]). In the thin capacitor model h/w — > 0, 
thus the force density vanishes. Thus, the total force acting on the wall is zero, and the 
above paradox is resolved. 

Concluding this Section we would like to discuss the applicability of Eq. flH]) to the case of 
a partially compensated zig-zag wall. To calculate the force acting on such wall, one should 
know the electric field and ferroelectric part of polarization on both sides of the wall and 
apply Eq. (jSJ) locally. In general, this is a tough mathematical problem. However, a simple 



w 

FIG. 2. Schematic of the zig-zag domain wall in a thin ferroelectric capacitor. The polarizations 
inside the domains are perpendicular to the electrodes. The bottom electrode is grounded and 
the constant potential, —U, is applied to the top electrode. The electric field Eq is uniform. The 
resultant force acting on the zig-zag wall depends only on the coordinates zq and z n of the side 
points. This force represents the edge-effect, which can be neglected in the thin capacitor model. 
Neglecting the edge-effects, the resultant force acting on the domain wall is zero. 

answer can be obtained in the case where the zig-zag amplitude d is much smaller than the 
domain size L (see Fig. [3]). We can consider such a zig-zag wall together with a layer of 
thickness L w , d <C L w <C L containing it. This layer is equivalent to the flat wall: they 
have the same jump of polarization and the same free charge, which create the same electric 
field in the bulk of the domains. The virtual displacement of the walls will lead to the same 
changes of the charges on electrodes and the same changes in the fraction of the domains. 
Thus the total forces acting on the zig-zag wall and the flat wall should be the same in the 
limit of d/L <C 1. This enables us to use Eq. ([5]) for the calculation the force acting in the 
walls of this type. 

III. FORCE ON A DOMAIN WALL BETWEEN FINITE DOMAINS 

In the previous section we obtained the force acting on the domain wall, carrying an 
arbitrary free charge density erf. The calculation of such force requires the knowledge of this 
charge. An essential interesting feature of this problem is that it is dependent, in a unique 
way, on other parameters of the system. This makes it possible to calculate the force, acting 
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FIG. 3. Schematic of zig-zag domain walls in a ferroelectric capacitor. The polarizations inside 
domains are perpendicular to electrodes. The zig-zag amplitude d is much smaller than the domain 
size. This wall can be considered together with a layer around it of thickness L w (d <C L w <C L), 
which is shown with the dark region on the figure. The polarizations inside the bulk of the domains, 
the free electric charges on the walls, and the force acting on the walls are the same for the flat 
and zig-zag walls. 




FIG. 4. Schematic of the periodic structure of " head-to-head" and "tail-to-tail" 180-degree domain 
walls with period L. Domain walls are perpendicular to the polarization inside the domains P, 
which is different from the spontaneous polarization due to the depolarizing field. 



on a charged wall in a domain pattern, due to the application of an external electric field. 
We will illustrate this for the case of a periodic domain pattern with charged walls (see Fig. 
Hj) in a uniaxial ferroelectric with the second-order phase transition. We will be interested 
in the situation, typical for ferroelectrics 17 , where the screening of the bound charge in the 
wall is in the nonlinear regime, i.e. the carrier concentration inside the domain wall is much 



higher, than the homogeneous carrier concentration in the material. We will consider a 
pattern with domains with the equal thicknesses x = L/2, where L is the structure period, 
having planar walls perpendicular to the spontaneous polarization. 

As was shown in Refs.— ,— , the charged wall in a finite sample is not completely com- 
pensated. The net charge on the wall is proportional (due to the Gauss law) to the electric 
field inside the adjacent domains. In the absence of external bias field, in a pattern con- 
taining charged domain walls, this field is not zero. For the considered situation, where 
the screening regime in the wall is nonlinear this field can be readily estimated. Indeed, in 
this regime, the free carrier density, needed for the polarization screening, can be provided 
exclusively by a band bending. Obviously, such band bending implies a voltage difference 
between the "head-to-head" and "tail-to -tail" domain walls, which is about E g /q, where 
E g is the band-gap of the ferroelectric, q is the absolute value of the electron chargeii^. 
Finally, the corresponding electric field inside a domain of the thickness L/2 is 

< 16 > 

It can actually also be viewed as a partially compensated depolarizing field. Due to this 
field the polarization inside the domains is reduced, compared to the nominal spontaneous 
polarization of the ferroelectric. For a ferroelectric with second-order phase transition, using 
equation of state E = aP + (3P 3 , the reduced polarization P can be found as the negative 
solution to the following equation: 

a p + f3pz = 2 -^. (17) 
qL 

A. Pressure on the wall between domains of equal sizes 

For the configuration considered, the force component parallel to the wall, given by Eq. 
(JSJ), is zero, whereas the pressure, given by Eq. (JTj), is non-zero. Considering domain wall 
between domains 1 and 2, one can rewrite Eq. © as follows: 

p = $ 2 _$ 1 _I( jBl + jB2 )(P 2 _P 1 _ (Tf ) ) (i 8 ) 

where $ = (a/2)P 2 + (/3/4)P 4 . Here, the electric field and polarization are considered 
positive if corresponding vectors are directed along the z axis (Fig. Hj), i.e. from domain 1 to 
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domain 2. It is clear from the structure of Eg. (1181) that the absolute values of the pressure, 
acting on the head-to-head and tail-to-tail walls, are equal, whereas the corresponding forces 
are acting in opposite directions. Thus, it suffices to calculate the pressure on the head-to- 
head walls between domains. 

The Poisson equation for such configuration has the form (see Eq.flU) for the definition of 
the electric displacement) : 

a i = e h E 2 + P 2 -e h E 1 -P 1 , (19) 

where electric fields inside the domains can be calculated using the equations of state: 

E, = aPi + (3P? } (20) 

E 2 = aP 2 + f3Pl (21) 
The external electric field E defined as E = U/h can be expressed in terms of E 1 and E 2 : 

E = {E 1 + E 2 )/2. (22) 

Using Eqs. ( |T9l) -(l22l) we can present the polarizations inside the domains as a Taylor 
expansion with respect to Eq as follows: 

P^P + xEo + vEl (23) 

P 2 = -P + X Eo-vEl (24) 



where 



X = (25) 

a + 3f3P 2 



is the susceptibility of the ferroelectric at P = P and 

= Z 3/3P = -3/3P X 2 

V (a + 3pP 2 ) 2 (a + 3/3P 2 + l/e h ) l/x + l/e b ' 

Substituting polarization from Eqs. (|23|) . ( 124|) into Eq. (fl8|) and keeping terms up to third 
power with respect to E , one finds: 

p = £l(2P + <r f )E - 2f3P X 3 E 3 , (27) 

where 

e f = e h + x (28) 
11 



is the permittivity of the ferroelectric at P = P. 
From the Poisson equation, Eq. ( fl9l) . we can find: 



2P + o-f = 2e h E, (29) 

where E is the electric field inside the domain, which in turn can be found from Eq. (fl6|) . 
Thus we obtain: 

2P + a f = ^. (30) 
qL 

Using Eq.( l30|) . we rewrite Eq. ( |27|) as follows: 

V = - 2(3P X 3 El (31) 

Note that in Eq. ( 13T]) the x an d £ f should be calculated for the reduced values of the 
polarization inside domains. 

It is instructive to compare the part of the pressure, linear in E , given by Eq.f l3ip for the 
case where the polarization inside the domains is close to the bulk spontaneous polarization 
Pq with the pressure on the neutral domain wall 2PqEq. Physically, the condition of the 
polarization close to Pq means that the depolarizing field inside domains, E, is much smaller 
than the thermodynamic coercive field E coer = (2/(3V3))\a\P . In this case, neglecting the 
difference between P and Pq and using Eqs. ( 1311) . (1251) one finds p/(2P E ) « 0.2E/E coer . 
Since, in the situation considered, the depolarizing field inside the domain E is much smaller 
than -E coer , the pressure on the charged wall is expected to be much smaller than the pressure 
on the neutral wall. 

The situation, where the depolarizing field approaches the thermodynamic coercive field, 
corresponds to a polarization instability inside domains manifesting itself in a divergency of 
the permittivity £f. This formally implies (via Eqs. (!3~T!) and (flET) ) an unlimited increase of 
the pressure acting on the wall when the domain period L approaches a certain critical value. 
This academically interesting situation is, however, hardly reachable experimentally^. For 
this reason, in this paper we will further discuss only the situation where the ferroelectric 
inside domains is far from the aforementioned instability. On the practical level, this means 
that in further estimates we can neglect the difference between the value of the polarization 
inside domain, P, and the spontaneous polarization of the ferroelectric Pq. 

12 



B. Polling a sample with charged domain walls 



In this section we will consider the problem of poling a sample, where the periodical 
structure with charged domain walls shown in Fig. H]was formed. We will consider the case 
of a small external field, where we can consider only the linear in Eq part of the pressure. 
Thus, keeping only the linear term in Eq. (l3Tj) . we find this pressure in the form: 

P = ^E a . (32) 
qL 

For the model developed in the previous section and used in the present one, the free 
carrier concentration in the domain wall is close to metallic. This means that, for a realistic 
concentration of the trapping centers, the main part of the free carriers in the wall is not 
trapped, having high mobility and, thus, being able to easily follow the domain wall in its 
motion. A realistic situation, which we consider, is that where the time needed for the 
free-charge to transfer between the domain walls (which depends on the conductivity of the 
material and domain size) is large in comparison to the time of poling. That means that in 
our model, when a domain wall moves away from its position in a periodic pattern, the free 
charge which its carries does not change. 

Defects, which always exist in real materials, lead to the domain wall pinning. The 
pinning acts as dry friction in mechanics. When we apply a small pressure on the domain 
wall, its position is stabilized because of the pinning. To start the domain wall motion, the 
applied pressure, given by Eq. (l32l) . should be equal to the maximal "stabilizing" pressure of 
the pinning p pinning : 

P Ppinning' 

(33) 

Equations ( 133]) and ( 1321) enable us to find the electric field needed to start the motion of 
the walls from the initial position in the equidistant domain pattern. To start the domain 
wall motion, it is necessary, but not sufficient for poling. To be sure that condition ( 133|) 
corresponds to the onset of successful poling, we should check if the pressure p is greater or 
equal to the pinning pressure at any position of the wall after it starts moving. If the pining 
pressure is independent of the position of the domain wall, successful poling will occur if, 
at given external electric field, the pressure acting at the wall in a shifted position is larger 
than that in the original position. In reality, one expects the pining pressure to be larger in 
the original position because of the defect accumulations at the wall during the time before 
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the poling; the situation observed experimentally 7 . It is clear that for this situation the 
aforementioned criterium of successful poling holds as well. 

Let us check if the criterium of successful poling is met for our system, shown in Fig. HI 
Since the absolute value of the pressure acting on the head-to-head and tail-to-tail walls are 
equal, whereas the corresponding forces are acting in the opposite directions, the period of 
the pattern will not be affected by the application of an electric field. Thus, when discussing 
poling we can address just one head-to-head wall moving between two tail-to-tail walls 
spaced by a time- independent distance L (see Fig. 0J. 

Let us first consider how the polarizations P±, P 2 and electric fields Ei, E 2 inside the 
domains will change after a domain wall displacement from initial position. Differentiating 
the Poisson equation, Eq. f[T9"j) . and taking into account that the free charge on the wall 
remains constant we get: 

dE 2 \dP 2 ( 8E X \dP 1 



(35) 



dP 2 J dx \ dPi J dx 
Here dE/dP « x' 1 ~ % \ thus e h {dE/dP) « e b /e f < 1 and we can conclude that 

dP 1 dP 2 

dx dx 

The sign of the derivatives in Eq. (I35j) can be obtained from the relation between average 
field Eq applied to the ferroelectric and fields inside the domains E\ and E 2 . 

E x x + E 2 (L -x) = E L. (36) 

Differentiating this equation we obtain: 

+ -.)--* + ft. (37) 

In the case of head-to-head wall considered here (see Fig. Hj) E x < 0, E 2 > 0, implying 
that the expression in the right hand side in Eq. (1371) is positive. Thus, the left hand side 
expression in Eq. f )37j) is also positive and we obtain 



^^^-^-dP^^dP,^^-^ >0 ' (38) 
Here dE\jdP\ and 8E 2 /dP 2 are inverse susceptibilities (which are positive); domain widths 
x, L — x are also positive. Thus, from Eq. (|38|) using Eq. (|35|) it follows 

dPx/dx w dP 2 /dx > 0. (39) 
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Below we will consider the positive displacement of the domain wall, i.e. the size of the 
domain 1 x is increasing during this displacement. This situation corresponds to the external 
field applied from domain 1 to domain 2 (E > 0). For this case, using Eq. (j3"91) we can 
conclude that polarization in domain 1, Pi > 0, and it will increase (the absolute value 
of depolarizing field E\ in this domain will decrease) while in the domain 2, P 2 < 0, and 
its absolute value will decrease (the absolute value of the field E 2 will increase). We can 
summarize the result for polarization change in a compact form as: 

\P 2 \ <P< <P , (40) 

| £i| < |E| < |E 2 |. (41) 
Now, using Eqs. (fl8l) . (fl9l) . we can obtain the pressure on the domain wall, when it is in 



p = $ 2 -$ 1 + |(E 2 -P 2 ). (42) 



an arbitrary position, as: 



In view of Eq. ( 14T1) when the domain wall moves in the positive sense of the z-axis (x > 1/2), 
E 2 increases, E\ decreases, thus the last term in Eq. (14*21) increases. Thus, to show that 
the pressure increases when the wall is displaced from the central position it is enough to 
check that term $ 2 — increases. Using Taylor expansion of polarization similar to Eqs. 
(!23|) - (l25|) and keeping only linear terms with respect to E we find: 

$ 2 - $ x = G(P U P 2 ) + if (A, i S 2)E , (43) 

where 

G(Pi, P 2 ) = \/3{P? ~ i 3 |)(2P 2 - Pi - A 2 ), (44) 



H(P lt P 2 ) = " 2l f 2 ~'° ; - " U f 1 ~ ro) . (45) 



P 2 (Pj - PI) Px(P? - Pp 2 
3P| - P 2 3P 2 - P 2 
Here Pi, P 2 , Ei, P 2 are polarizations and electric fields in non-symmetric situation (x ^ 
L — x) at zero external field in the domains 1 and 2, respectively. 

In the case of the zero external field Eq. fj4*0|) leads to P 2 < P\ < Pq, and it immediately 
follows that G(Pi,P 2 ) > 0. Thus if we displace the domain wall from the central position, 
and the external field is zero, the force will act on the wall, which is parallel to the displace- 
ment. Thus the equidistant periodic structure with charged domain walls is unstable. It 
corresponds to a maximum of the energy. In real systems, the position of charged domain 
walls are stabilized because of the pinning on defects, discussed above. 

15 



Now we will analyze the second term in Eq. ( 143]) . which is proportional to Eq. Using 
Eqs. (@5D, flMD we find 



dH „ 6(P 2 P 2, Pq(p? + Pi - P 2 ) + ^P?Pj dp 

dx ~ 1 1 2 J (3P| - P 2 ) 2 (3P 2 - P 2 ) 2 9k ' { ' 

From Eq. (146]) using Eqs. (139]) . (f4T)]) and taking into account that Pi ~ P , P 2 ~ P we 
have dH/dx > 0. 

Thus after we displaced domain wall from the central position by application of the 
external electric field, the pressure acting on the wall will increase, and the criterium of 
successful poling is met. This means that the switching is governed by the condition in the 
original position of the domain wall given by Eqs. ( 133]) and ( 132]) . 



C. Nonlinear part of the pressure 

The nonlinear part of the pressure Eq. (1271) is directed in the opposite direction to the 
linear part, and, for the structure with large domains where P — > Pq, it leads to the nonlinear 
pressure obtained as^: 

p=-(x 2 /Po)El (47) 

If the domain wall is well compensated and the applied field is high enough, the nonlinear 
pressure can play the decisive role and resultant pressure can act opposite to the linear 
one. Thus, if we abruptly apply high voltage, the domain wall may move in the direction, 
opposite to the direction of motion at low voltage. Let's estimate the electric field at which 
the values of the nonlinear and linear pressure are equal. Setting P « Pq from Eq. ( }3T|) one 
can find: 

E = \J 4P g P / qLx- (48) 

For the typical values E g = 3 eV, P = 0.1 C/m 2 , \ = 10 3 x £ o> L = 10 A™ (this is about 
the period of the structure observed by Shur-) one can check that the corresponding field 
inside the domain is far from the thermodynamic coercive field. Thus Eq. (148]) can be used 
and it leads to Po ~ 40 kV/cm. This field is decreasing with increasing the domain size. 
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IV. DISCUSSION AND CONCLUSIONS 



We obtained a general formula for the local force acting on the domain wall with arbitrary 
free charge on it. In general, this force in not perpendicular to the domain wall. The force 
acting on the domain wall depends on the compensating free charge on it. The better the 
compensation, the smaller the force acting on the domain wall. The force acting on the fully 
compensated wall is zero. 

The compensating free charge on the domain wall depends on the characteristic size of 
the domain structure. We considered the periodic domain structure with parallel "head-to- 
head" and "tail-to-tail" walls, and obtained the force acting on the domain wall as a function 
of the structure period. Smaller domain size leads to a larger force. For the structure which 
is not very dense the force acting on the charged domain wall is much less than the force 
acting on the neutral domain wall. Thus poling of the sample with such a structure is 
difficult in comparison with poling of the sample with neutral domain walls. This result is 
in good agreement with the experimental data obtained by Shur et al.— 

We also showed that in an ideal crystal, the periodic structure with charged domain walls 
is unstable. It corresponds to the maximum of the energy. Nevertheless, in real materials it 
can be stabilized due to pinning provided by defects. 
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Appendix A: Derivation of the local force density acting on the element of the 
domain wall 

We analyze the local pressure on the ferroelectric domain wall in a system shown in Fig. 
EJ We consider that inside a ferroelectric material with ferroelectric part of polarization 
there exists a closed domain of the material in a different polarization state P^ ■ The 
closed boundary splits the ferroelectric into two domains, which have volumes and 
V( 2 \ respectively. Inside domains there are conductors, which carry charges and , 
respectively, which are at electric potentials tp^ and <f^i respectively. We consider that 
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FIG. 5. General system of a ferroelectric in contact with two conductors. Domain wall Sw separates 
the ferroelectric into two domains of volumes and V ^ . There are two conductors, which carry 
charges and and have electric potentials tp^' and <p& within each domain. The quantities 
within each domain are the ferroelectric part of polarization P^ 1 ' and JR , electric field and 
E^ and electric potential (p^ and ip( 2 \ Symbol 5u stands for the virtual displacement of the 
domain wall and symbol / stands for the local force density of external sources on the domain wall 
in order to keep the system in equilibrium. 

the bound charges due to discontinuous change of polarization at the domain wall Sw are 
partially compensated by free charges of surface density <jf. The charges on conductors and 
on the domain wall produce electric fields E? and within each domain. Symbols tp^ 
and ip^ stand for the electric potentials within each domain. 

In order to obtain a general formula for the local pressure on domain wall, we follow the 
principle of virtual displacements, formulated in SectJTTJ Instead of Eq. (jSJ) we will use the 
thermodynamic function G in a more general form: 

dV-X>«Q« (Al) 

i 

where the first term represents the energy associated with the ferroelectric part of polariza- 
tion, $(-P), and the electric field energy integrated over the volume V of ferroelectric, and 
the second term represents the subtracted work of the electric sources. 

In what follows it is convenient to transform the work of electric sources into volume 



G 



$(P) + -e h E 2 
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integrals: 

Q ( e = -V { e [ D<p ni dS=[ E^D^dV+f ^D^n t dS, (A2a) 
Jsg> Jvw Js w 

<P$Q$ = -<P$ f Df\dS=! E^D^dV-f ^D^n t dS, (A2b) 
Jsf Jvw Js w 

where and are the vectors of electric displacement in domains and V^ 2 \ 

respectively, and where we considered the absence of free charges, i.e. div-D, = 0, inside 
domains V^ 1 -* and and vanishing the surface integral over the domain at infinity. 
Orientations of the normal vectors rij are indicated in Fig. Since the bulk quantities, e.g. 
<pDi, are generally discontinuous at the domain wall Sw, it will be useful to denote their 
jump by, e.g. 

[[pD i ]]=^Df^-^Df\ (A3) 

If we employ the expression for electric displacement D i = P i J r s^Ei, the thermodynamic 
function G can be expressed in a form: 



G 



+ 



+ 



y(i) 



y(2) 



Sw 



'^) {P m ) _ l _ £hE f) E f)_ E f) P p 

P (P l + e h E i )]] mdS. 



dV + 
dV + 



(A4) 



Now our task is to express the variation of thermodynamic function 5G, which is produced 
by the virtual displacement of domain wall. When doing this, we have to keep in mind that 
the bulk quantities tp, Ei, and Pj are not independent variables. The relation between the 
electric field and the electric potential is given by Ei = —dip/dxi and one has to consider 
the continuity of electric displacement at the domain wall, i.e. [[Pi + e^Ei]] rii = Of. These 
additional conditions can be introduced into our treatment by adding new variables \ and 
H according to the method of Lagrange multipliers. Using this approach we obtain a new 
thermodynamic function L = L(ip, Ei, Pi, Aj, //), where all bulk variables in the both domains 
are considered independent: 



+ 



y(i) 



y(2) 



Sw 



- ^ h E^E^ - E^PV + A^ ( E? ] + 



(l)o(l) , \(1) 



,(1) 



dip 



(i) 



$ (2) (^ (2) ) - ^EfEf - E^pf* + Af ( EP 



dS. 



2 
1 
2 

[[if (Pi + e h Ei)]] m + fi ([[Pi + e h Ei]] m - a f ) 
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dxi 
dtpW 
dxi 



dV + (A5) 
dV + 



The virtual displacement Su of domain wall, which in general is not normal to the wall 
surface, shown in Fig. [5] produces the variation of all independent quantities Sip, SE iy SPi, 
5Xi, and Sfi; except the external force / on the domain wall and the electric potential on 
conductors (p^ and (p^\ which are constant in the system. The variation of quantities 
results in the variation of thermodynamic function SL. In equilibrium the variation SL 
equals the work produced by external force during the virtual displacement of domain wall 

SW = -f fiSuidS, (A6) 

where is the local force density per unit area of the domain wall, S\y is the area of domain 
wall. 

In order to express the variation SL, it is useful to consider following formula for the 
variation of volume integral 

S\ [ gdv\= [ SgdV± [ g5u k n k dS, (A7a) 
{Jv J Jv Js w 

where g is an arbitrary bulk quantity and plus and minus correspond to volume and 
V^ 2 \ respectively. The physical interpretation of the above formula is that the variation 
of volume integral includes the contribution due to variation of bulk quantity Sg and the 
contribution due to the volume change of domains produced by the virtual displacement of 
domain wall Su. Similar formula can be written for the variation of surface integral over the 
domain wall 

S\ f gds\= [ SgdS+ [ p-Su k dS, (A7b) 

lJs w ) Js w JSw ® x k 

where the first term on the right-hand side represents the contribution due to variation of 
bulk quantity Sg on the domain wall and the second term represents the contribution due 
to change of the domain wall position during the virtual displacement. 
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Applying Eqs. ( 1A7j) to function L given by Eq. (1A5|) . we obtain 



5L 



5P 



(i) 



-P 



(i) 



ap 



5£> 



(i) 



+ 



y(2) 



6P 



(2) 



^ A ( 2) + 5X m 



E 



P^+e^-X^ 
dV+ 



(A8) 



dxi 



SE} 



(2) 



p (2) + £b p ? (2) - a! 2) 



+ 



+ 



Sw 



5ip 



(2) 



2) 



(i) 



+5/i([[P + £ b P i ]]n i + 



9xj 



$(P) - -StEiEi - E^ + Xi [Ei + 



dip 



([[<p (Pi + e h Ei) m)} + n ([[Pi + e h E^] m - ay)) J \ dS. 



The above formula can be further transformed using following integral expressions: 



I (1) / dS^_ \ dy+ f fd^\ {l)dv = I (1) 
Jvw \ dxi J J vW y dxi J J Sw 



where we consider that, during the virtual displacement of domain wall, the electric potential 
on conductors is kept constant by the electric sources, i.e. the variations of electric potentials 
on conductors are zero, and that the surface integral over domain is vanishing at infinity. 
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Applying Eqs. ( 1A9j) to formula (1A8j) . we readily obtain the variation SL in a form 



y(i) 



6P 



(i) 



OP 



E 



(i) 



5P 



(i) 



d\ (1) 



+ 5A 



(i) 



dV+ 



(A10) 



+ 



V(2) 



<5P 



(2) 



E 



(2) 



5$(2) 



rj-. (2) 

ax* 



,(2) 



5E\ 

ET' : 



<9xi 

p^ + s^r-A; 



(2) , ^ (2) 



(2) 



P^+e^-A 



+8fi([[P i + s h E i ]]n i -a f ) + 



-5u k 

_d_ 

dxi 



$(P) - ~£ b P*Pi - PiP + A 4 ( P + 



([[<? (P + e h E { ) m]] + /i ([[P + e h Ei}} m - a f ))j J dS, 



which can be combined with the work of external pressure SW during virtual displacement of 
domain wall given by Eq. (1A6|) . Employing the principle of virtual displacements 5L = 5W 
yields the formula for the local force density of external sources on the domain wall 



fk 



$(P) - -e^Ei - E^ + \ i [E i + 



d 



dip 
dxi 



(All) 



dxi 



{[[if (P + e h Ei) m}} + fi ([[Pi + e h Ei}} m - a f )} , 



equations of motion 



<9$ 
OP 



(A12) 



electric displacement Aj = e^Ei + Pi, Gauss' law for electric displacement dXi/dxi = 0, 
relationship between electric field and electric potential Pj = —dip/dxt, continuity of electric 
potential at the domain wall fi = ip^ = ip^ 2 \ and continuity of electric displacement at the 
domain wall [[Pi + £bP«]] n % — °> Combining the above expressions, the formula for the 
local force density of external sources on the domain wall can be written in a form 



fk 



$(P) - -e h EiEi - E t P t 



n k + [[E k (Pi + e h Ei)\] 



Hi. 



(A13) 
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Considering the continuity of tangential components of electric field at the domain wall 
[[E^j]] = \[Ei — (EkUkjUi]] = and using the algebraic identity 

[[fg]] = f[[g]] + [[/]]?, (am) 

where / = (/W 4- fi 2 )^ ji j s the average of bulk quantity at the opposite sides of the domain 
wall, the general formula for the local force density of external sources on the domain wall 
can be further simplified to the form: 

f k = [MP)]) n k - % [[Pi]} n k + E k a t . (A15) 
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